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A search for diphoton events with large missing transverse energy is presented. The data were 
collected with the ATLAS detector in proton-proton collisions at y/s = 7 TeV at the CERN Large 
Hadron Collider and correspond to an integrated luminosity of 3.1 pb . No excess of such events 
is observed above the Standard Model background prediction. In the context of a specific model 
with one Universal Extra Dimension with compactification radius R and gravity-induced decays, 
values of 1/R < 729 GeV are excluded at 95% CL, providing the most sensitive limit on this model 
to date. 



PACS numbers: 14.80.Rt, 13.85.Rm 

In the Standard Model (SM), the production in proton- 
proton (pp) collisions of diphoton (77) events with large 
missing transverse energy (Kj? lss ) is mainly due to W/Z+ 
77 processes. Taking into account the branching ratios 
of W/Z decays including at least one neutrino, the cross 
sections are only a few femtobarns for 7 TeV pp collisions. 
In contrast, some new physics models predict much larger 
77 + E^ iss rates. This Letter reports the first 77 + E^ iss 
search with LHC data, using data recorded with the AT- 
LAS detector. The results are interpreted in the context 
of a Universal Extra Dimension (UED) model. 

UED models [l| postulate the existence of additional 
spatial dimensions in which all SM particles can prop- 
agate, leading to the existence for each SM particle of 
a series of excitations, known as a Kaluza-Klein (KK) 
tower. This analysis considers the case of a single TeV -1 - 
sized UED, with compactification radius R. The masses 
of the states of successive levels in the tower are sepa- 
rated by ~ 1/R. For a given KK level, the approximate 
mass degeneracy of the KK excitations is broken by ra- 
diative corrections 0. The lightest KK particle (LKP) 
is the KK photon of the first level, denoted 7*. At the 
LHC, the main UED process would be production via 
the strong interaction of a pair of first-level KK quarks 
and/or gluons Q , which would decay via cascades involv- 
ing other KK particles until reaching the LKP at the end 
of the decay chain. If the UED model is embedded in 
a larger space with N additional eV _1 -sized dimensions 
accessible only to gravity 0| , the LKP could decay grav- 
itationally via 7* — >• 7 + G Q, where G represents one 
of a tower of eV-spaced graviton states. With two decay 
chains per event, the final state would be r yy + E^ lss + X, 
where E™ lss results from the escaping gravitons and X 
represents SM particles emitted in the cascade decays. 

The UED model considered is defined by specifying R 
and A, the ultraviolet cut-off used in the calculation of 
radiative corrections to the KK masses. This analysis 
treats R as a free parameter and, following the theory 
calculations @, sets A such that AR = 20. For 1/R = 
700 GeV, the masses of the first-level KK photon, quark 



and gluon are 700, 815 and 865 GeV, respectively Q. 
The 7* mass is insensitive to A, while other KK masses 
change by typically a few percent when varying AR in 
the range 10 — 30. The gravitational decay widths of the 
KK particles are set by N and Md, the Planck scale in 
the (4 + A)-dimensional theory. For the chosen values of 
N = 6 and M D = 5 TeV, and provided 1/R < 1 TeV, 
the LKP is the only KK particle to have an appreciable 
rate of gravitational decay. The same parameter values 
were used in the only previous study of this model, in 
which the DO experiment excluded at 95% CL values of 
l/i?< 477 GeV 0. 

Monte Carlo (MC) signal samples were produced for 
a range of 1/R values using the implementation [(| of 
the UED model in PYTHIA || version 6.421, and using 
the MC09 parameter tune [s| . The MC samples were pro- 
cessed through the ATLAS detector simulation [lOj based 
on GEANT4 llf . In addition to the two high transverse 
energy (Et) photons and large E™ lss , the signal events 
typically include several high-£>r jets due to the cascade 
decays, with the Et spectrum of the leading jet peaking 
at w 100 GeV for 1/R = 700 GeV. 

The ATLAS detector [l2j is a multipurpose particle 
physics apparatus with a forward-backward symmetric 
cylindrical geometry and nearly 47r solid angle cover- 
age. ATLAS uses a Cartesian right-handed coordinate 
system, with the nominal collision point at the origin. 
The anticlockwise beam direction defines the positive 
z-axis, while the positive x-axis points from the colli- 
sion point to the centre of the LHC ring and the posi- 
tive y-axis points upward. The angles <fi and 9 are the 
azimuthal and polar angles. The pseudorapidity is de- 
fined as 77 = — In [tan (0/2)}. Closest to the beamline are 
tracking detectors which use layers of silicon-based and 
straw-tube detectors, located inside a thin superconduct- 
ing solenoid that provides a 2 T magnetic field, to mea- 
sure the trajectories of charged particles. The solenoid is 
surrounded by a hermetic calorimeter system. A liquid- 
argon (LAr) sampling calorimeter is divided into a central 
barrel calorimeter and two end-cap calorimeters, each 
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housed in a separate cryostat. Fine-grained LAr electro- 
magnetic (EM) calorimeters, with excellent energy res- 
olution, provide coverage for \rj\ < 3.2. In the region 
1 77 1 < 2.5, the EM calorimeters are segmented into three 
longitudinal layers and the second layer, in which most of 
the EM shower energy is deposited, is divided into cells 
of granularity of A77 x A</> = 0.025 x 0.025. A presam- 
pler, covering \r}\ < 1.8, is used to correct for energy 
lost upstream of the calorimeter. An iron-scintillator 
tile calorimeter provides hadronic coverage in the range 
1 77 1 < 1.7. In the end-caps (\r]\ > 1.5), LAr hadronic 
calorimeters match the outer |ry| limits of the end-cap 
EM calorimeters. LAr forward calorimeters provide both 
EM and hadronic energy measurements, and extend the 
coverage to |?7| < 4.9. Outside the calorimeters is an 
extensive muon system including large superconducting 
toroidal magnets. 

The reconstruction of photons is described in detail in 
Ref. [3]. To select photon candidates, EM calorimeter 
clusters were required to pass several quality criteria and 
to lie outside problematic calorimeter regions. Photon 
candidates were required to have \r]\ < 1.81 and to be 
outside the transition region 1.37 < \rj\ < 1.52 between 
the barrel and the end-cap calorimeters. The analysis 
uses a "loose" photon selection, which includes cuts on 
the energy in the hadronic calorimeter as well as on vari- 
ables that require the transverse width of the shower, 
measured in the second EM calorimeter layer, be consis- 
tent with the narrow width expected for an EM shower. 
The loose selection provides a high photon efficiency with 
modest rejection against the background from jets. 

The reconstruction of E™ lss is based on topological 
calorimeter clusters 14| with |r7| < 4.5 that are seeded 



by any cell with energy higher than four times its noise 
level. In an iterative procedure, the cluster grows by in- 
cluding all neighbouring cells with energy higher than 
twice the noise, plus all cells neighbouring the bound- 
ary of this three-dimensional collection. Each cluster is 
classified as EM or hadronic, depending on its topology, 
and the cluster energy is calibrated to correct for the 
non-compensating calorimeter response, energy losses in 
dead material, and out-of-cluster energies. Events re- 
constructed with large E^ lss were studied in detail with 
early data 453- Rare background events with large trans- 
verse energies, unrelated to the collision and concentrated 
in a few cells, due mainly to discharges and noise, have 
been observed. Cuts were applied to eliminate such back- 
grounds, rejecting less than 0.05% of the selected events 
while having a negligible impact on the signal efficiency. 

The data sample was collected during stable beam pe- 
riods of 7 TeV pp collisions at the LHC, and corresponds 
to an integrated luminosity of 3.1 pb _1 . The events se- 
lected had to satisfy a trigger requiring at least one loose 
photon candidate with Et > 20 GeV, and had to con- 
tain at least one reconstructed primary vertex consistent 
with the average beam spot position and with at least 



three associated tracks. The trigger and vertex require- 
ments are ~ 99% efficient for signal MC events. The 
presence of multiple pp collisions within the same bunch 
crossing, known as "pileup" , can be analysed by examin- 
ing N vtx , the number of reconstructed primary vertices 
in each event. In this data sample, the average value 
of N vtx was w 2.1. The MC signal samples included the 
simulation of pileup and were weighted to match the N vtx 
distribution observed in data. 

Events were retained if they had at least two photon 
candidates, each with Et > 25 GeV. In addition, a pho- 
ton isolation cut was applied, wherein the Et in a radius 
of 0.2 in the rj-cj) space around the center of the cluster, 
excluding the cells belonging to the cluster in a region 
corresponding to 5 x 7 cells in 77 x in the second layer 
of the EM calorimeter, had to be less than 35 GeV. This 
requirement had a signal efficiency greater than 95% but 
rejected some of the background from multijet events. 
An event in which each of the two photon candidates 
satisfied the loose photon cuts was considered a 77 can- 
didate event. An independent "misidentified jet" control 
sample, enriched in events with jets misidentified as pho- 
tons, was defined as those events where at least one of the 
photon candidates did not pass the loose photon identi- 
fication. After all cuts, the 77 and misidentified jet sam- 
ples totaled 520 and 7323 events, respectively. Figure Q] 
shows the Et spectrum of the leading photon for the 77 
candidates and for UED l/R = 700 GeV MC events; the 
UED spectrum extends to much higher Et values. 
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FIG. 1. Et spectrum of the leading photon for the 77 candi- 
date sample and for UED 1/R = 700 GeV MC events (nor- 
malized to 100 times the leading order (LO) cross section). 

The background was evaluated entirely using data. 
Non-collision backgrounds, such as cosmic rays and 
beam-halo events, are reduced to a negligible level by the 
selection cuts. The main background source, referred to 
hereafter as QCD background, arises from a mixture of 
SM processes including 77 production, and 7 + jet and 
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multijet events with at least one jet misidentified as a 
photon. With the loose photon identification, it is ex- 
pected that 7 + jet and multijet events dominate, with 
only a small 77 contribution. The misidentified jet sam- 
ple provided a model of the E™ iss response for events 
with jets faking photons. The response for 77 events was 
modeled using the E™ lss spectrum measured in a high 
purity sample of Z — >• ee events, selected by a combina- 
tion of kinematic cuts and electron identification require- 
ments The £"™ lss spectrum for Z — > ee events, which 
is dominated by the calorimeter response to two genuine 
EM objects, was verified in MC to model the i?™ lss re- 
sponse in SM 77 processes, despite their kinematic differ- 
ences. As shown in Fig. [21 Z — > ee events typically have 
somewhat lower Kj? lss values than events of the misiden- 
tified jet sample, as expected since the presence of jets 
faking photons should result in a broader E™ lss distri- 
bution. The spectrum for the 77 candidates, which for 
low E™ aa is dominated by the QCD background with 
an unknown mixture of events with zero, one, and two 
fake photons, lies between these two samples. The £"™ ISS 
spectrum of the total QCD background was modeled by a 
weighted sum of the spectra of the Z — > ee and misiden- 
tified jet samples. The QCD background was normalized 
to have the same number of events as the 77 candidate 
sample in the region E™ laa < 20 GeV, where any UED 
signal contribution can be neglected. The relative contri- 
butions of the Z — > ee and misidentified jet samples were 
determined by fitting the QCD background shape to the 
spectrum of the 77 candidates in this same low 
The fraction attributed to 77 production, 
as modeled with the Z — ¥ ee distribution, was determined 
to be (36±22)%. The search result is not very sensitive to 
the exact composition of the QCD background, and the 
fit error was used to determine systematic uncertainties 
on the background prediction. 

A small additional background results from W — )• ev 
events, which have genuine E™ laa and which can pass the 
selection if the electron is misidentified as a photon and 
the second photon is either a real photon in W-f events 
or a jet faking a photon in W + jets events. A high 
purity sample of inclusive W — > ev events was selected 
by a combination of kinematic and electron identifica- 
tion cuts [3] ■ Requiring in addition a loose photon with 
> 25 GeV, a "W + 7" sample of only 5 events was se- 
lected. Accounting for the probability for an electron to 
be misidentified as a loose photon, as determined using 
the Z — > ee sample, the total background contribution 
due to W — > ev events was then estimated to be only 
w 0.4 events. Since the number of Wj events was too 
small to measure their E™ laa spectrum, a sample of W + 
jets events was used instead, requirin g a j et reconstructed 
with an anti-fcx clustering algorithm [161 ] with radius pa- 
rameter 0.4 and E{ > 25 GeV. The W(-> ev) + jets/7 
background contribution was then estimated by normal- 
izing the W + jets _E™ ISS spectrum to the expected total 
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FIG. 2. _B™ 1SS spectra for the 77 candidates, for the Z — > ee 
and misidentified jet samples used to model the QCD back- 
ground (each normalized to the number of 77 candidates with 
E^ iss < 20 GeV), and for the W(-> ev) + jets/7 background 
(normalized to its expected total of « 0.4 events). Variable 
sized bins are used, and the vertical error bars and shaded 
bands show the statistical errors. 



of « 0.4 events, as shown on Fig. 

Figure [3] shows the E™ lss spectrum of the 77 candi- 
dates, superimposed on the total background prediction, 
as well as example UED signals. Table Q] summarizes 
the number of observed 77 candidates, as well as the 
expected backgrounds and example UED signal contri- 
butions, in several E™ laa ranges. The QCD background 
dominates, and falls steeply with rising E™ laa , while the 
W — >■ ev background is very small, and flatter as a func- 
The UED signals would peak at large val- 
There is good agreement between the data 
and predicted background over the entire E™ laa range, 
with no indication of an excess at high E™ lss values. 

The signal search region was chosen to be E™ ss > 
75 GeV, before looking at the data, to obtain the best 
sensitivity to the UED signal. In the signal region, there 
are zero observed events, compared to an expectation of 
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0.32±0.16(sfat.)t u 'io( s y s ' : -) background events. The sys- 
tematic uncertainty was derived by studying variations of 
the background determination, including varying within 
its error the 77 fraction determined in the fit of the QCD 
background, varying the definition of the misidentified jet 
sample, and eliminating the photon isolation cut. 

The UED signal efficiency, determined from MC, in- 
creases smoothly from w 43% for 1/R = 500 GeV to 
w 48% for 1 /R = 700 GeV, with the lower efficiencies for 
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TABLE I. The number of observed 77 candidates, as well as the SM backgrounds estimated from data and expected UED 
signal for 1/R values of 500 and 700 GeV, given in various E™' BB ranges. The uncertainties are statistical only. The first row, 
for E™' aa < 20 GeV, is the control region used to normalize the QCD background to the number of observed 77 candidates. 



E^ iaa range 


Data 




Predicted background events 


Expected UED 


signal events 


(GeV) 


events 


Total 


QCD 


W{— > ev) + jets/7 


1/R = 500 GeV 


1/R = 700 GeV 


- 20 


465 


465.0 ±9.1 


465.0 ±9.1 




0.28 ±0.06 


0.02 ±0.01 


20 - 30 


45 


40.5 ± 2.2 


40.41 ± 2.17 


0.11 ±0.07 


0.45 ± 0.07 


0.03 ± 0.01 


30 - 50 


9 


10.3 ±1.3 


10.13 ±1.30 


0.16 ±0.10 


1.60 ±0.12 


0.08 ± 0.01 


50 - 75 


1 


0.93 ±0.23 


0.85 ±0.23 


0.08 ±0.05 


2.84 ±0.16 


0.14 ±0.01 


> 75 





0.32 ±0.16 


0.28 ±0.15 


0.04 ±0.03 


40.45 ±0.62 


4.21 ± 0.06 
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effects on the E!f lss reconstruction 14 1, including pileup, 
varying the cluster energies within the current uncertain- 
ties, and varying the expected 7?™ lss resolution between 
the measured performance and MC expectations, com- 
bine to give a 1% uncertainty on the signal efficiency. 
Finally, the 1% statistical error on the signal efficiency 
as determined by MC is treated as a systematic uncer- 
tainty on the result. Adding in quadrature, the total 
systematic uncertainty on the signal yield is 12%. 



TABLE II. Relative systematic uncertainties on the expected 
UED signal yield. For more details, see the text. 



Source of uncertainty 


Uncertainty 


Integrated luminosity 


11% 


Photon reconstruction and identification 


4% 


Effect of pileup 


2% 


E™ laa reconstruction and scale 


1% 


Signal MC statistics 


1% 


Total 


12% 



FIG. 3. E™ ma spectrum for the 77 candidates, compared 
to the total SM background as estimated from data. Also 
shown are the expected UED signals for 1/7? = 500 GeV and 
700 GeV. Variable sized bins are used, and the vertical error 
bars and shaded bands show the statistical errors. 



smaller 1/R due mostly to the 75™ lss > 75 GeV definition 
of the signal region. The various relative systematic un- 
certainties on the extraction of the UED signal cross sec- 
tion are summarized in Table UH including the dominant 
11% uncertainty on the integrated luminosity 17]. Un- 
certainties on the efficiency for reconstructing and identi- 
fying the 77 pair arise mainly due to differences between 
MC and data in the distributions of the photon identifi- 
cation variables, the need to extrapolate to the higher Et 
values (see Fig. Q]) typical of the UED photons, the im- 
pact of the photon quality cuts, varying the scale of the 
photon 7?t cut, and uncertainties in the detailed mate- 
rial composition of the detector. Together these provide 
a systematic uncertainty of 4%. The influence of pileup, 
evaluated by comparing MC samples with and without 
pileup, gives a systematic uncertainty of 2%. Systematic 



Given the good agreement between the measured £"™ lss 
spectrum and the expected background, a limit was set 
on 1/R in the specific UED model considered here. A 
Bayesian approach was used to calculate a limit based 
on the number of observed and expected events with 
£™ 1SS > 75 GeV. A Poisson distribution was used as 
the likelihood function for the expected number of signal 
events, and a flat prior was used for the signal cross sec- 
tion. Log-normal priors were used for the various sources 
of uncertainty, which were treated as nuisance parame- 
ters. It was verified that the result is not very sensitive to 
the detailed form of the assumed priors. Figure 0] depicts 
the resulting 95% CL upper limit within the context of 
the UED model considered, together with the LO UED 
cross section as a function of 1/R. The LO cross section 
was used since higher order corrections have not been 
calculated for the UED model. An uncertainty on the 
signal cross section due to parton distribution functions 
(PDF) was determined by comparing the predictions us- 
ing MRST2007 [3 PDF's with those from the full set of 
error PDF's of CTEQ6.6 [3]. The resultant uncertainty, 
namely ±8% essentially independent of 1/R, is shown by 
the width of the theory curve band. The observed 95% 
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CL exclusion region is 1/R < 729 GeV. The result de- 
pends weakly on the systematic uncertainties, and would 
only increase to 732 GeV if they were neglected. Chang- 
ing the E™ lss cut to 60 or 90 GeV would change the limit 
by only a few GeV. A cross check using a higher purity 
77 sample, achieved by requiring that both photons pass 
tighter identification cuts that reject more of the back- 
ground from jets, produced a consistent result. 
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FIG. 4. 95% CL upper limits on the UED production cross 
section, and the LO theory cross section prediction, as a func- 
tion of 1/R. The shaded band shows the PDF uncertainty. 

In conclusion, a search for 77 events with large E™ lss , 
conducted using a 3.1 pb _1 sample of 7 TeV pp collisions 
recorded with the ATLAS detector at the LHC, found 
no evidence of an excess above the SM prediction. The 
results were used to set limits on a specific model with 
one UED and gravity-induced LKP decays, excluding at 
the 95% CL values of 1/R < 729 GeV, and significantly 
surpassing the only existing experimental limit f7fl on this 
model. 
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